
Calhoun 

iniQiuiic^iul Ar{hiv« of tilt Mil vdl Poii^roduiit School 


Calhoun: The NPS Institutional Archive 
□Space Repository 



Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


2019-05 

Quantification of Extreme Event Statistics in 
Ship Design 


Rathmore, Udit 

Monterey California. Naval Postgraduate School 


http://hdl.handle.net/10945/63183 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 

Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


htt p://w ww. n ps. e du/l ib ra ry 


Callwuo is the Naval Postgraduate School's public access digital repository for 
research mate rials and institutiional publicatkios created by the NPS community. 
Calhoun is named for Professor of Mathematics Guy K. Caftiouo, NPS's first 
appointed — and published — schoteily author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 Univefsity Circle 
Monterey, California USA 93943 






REPORT DOCUMENTATION PAGE 


Form Approved 
0MB No. 0704-0188 


The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instruotions, searching existing data 
sources, gathering and maintaining the data needed, and completing and reviewing the colleotion of information. Send comments regarding this burden estimate or any other 
aspect of this collection of information, including suggestions for reducing the burden, to Department of Defense, Washington Headquarters Services, Directorate for Information 
Operations and Reports (0704-0188), 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302. Respondents should be aware that notwithstanding any other 
provision of law, no person shall be subject to any penalty for failing to comply with a collection of information if it does not display a currently valid 0MB control number. 

PLEASE DO NOT RETURN YOUR FORM TO THE ABOVE ADDRESS. 


1. REPORT DATE fDD-MM-YYyyj 2. REPORT TYPE 

05/10/2019 Master's Thesis 


4. TITLE AND SUBTITLE 

Quantification of Extreme Event Statistics in Ship Design 


6. AUTHOR(S) 

Udit Rathore 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Massachusetts Institute of Technology 
77 Massachusetts Ave., 

Cambridge, MA 02139 


3. DATES COVERED (From - To) 

09/24/2018 -05/10/2019 


5a. CONTRACT NUMBER 

N/A 

5b. GRANT NUMBER 

N/A 

5c. PROGRAM ELEMENT NUMBER 

N/A 


5d. PROJECT NUMBER 

N/A 

5e. TASK NUMBER 

N/A 

5f. WORK UNIT NUMBER 

N/A 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

N/A 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

N/A 


10. SPONSOR/MONITOR'S ACRONYM(S) 

N/A 


11. SPONSOR/MONITOR'S REPORT 
NUMBER(S) 

N/A 


12. DISTRIBUTION/AVAILABILITY STATEMENT 

DISTRiBUTiON A. Approved for pubiic reiease: distribution uniimited 


13. SUPPLEMENTARY NOTES 


14. ABSTRACT 

This thesis builds on previous work at the MIT Stochastic Analysis and Non-linear Dynamics (SAND) lab for the 
quantification of extreme events using wave groups. By separating the event probability from the physics models, we are 
able to capture rare events in ship motion and loading conditions for a modest computational cost. Improvements to the 
wave groups methodology ensured the slope and amplitude of the incident waves reflected the waves encountered in a 
given wave spectrum. The remaining discussion explores the value of a near-real-time risk analysis tools in reference to 
ship design and ship operations, with unique application to Navy and Commercial Vessels. 

15. SUBJECT TERMS 

Risk-based Ship Design, Wave Groups, Operational Efficiency, Extreme Event, Monte Carlo, Risk-related Design, 


16. SECURITY CLASSIFICATION OF: 

17. LIMITATION OF 

18. NUMBER 

19a. NAME OF RESPONSIBLE PERSON 

a. REPORT 

b. ABSTRACT 

c. THIS PAGE 

ABSTRACT 

OF 

PAGES 

MIT Academic Officer 






19b. TELEPHONE NUMBER (Include area code) 

u 

u 

U 

uu 

86 

617-253-4339 


standard Form 298 (Rev. 8/98) 
Prescribed by ANSI Std. Z39.18 















Quantification of Extreme Event Statistics 


in Ship Design 

by 

Uditbhan S. Rathore 
B.S., U.S. Naval Academy, 2012 

Submitted to the Department of Mechanical Engineering and 
System Design and Management Program in partial fulfillment of the requirements 

for the degrees of 

Master of Science in Engineering and Management 

and 


Naval Engineer 
at the 

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 

June 2019 

(c) Uditbhan S. Rathore, 2019. All rights reserved. 

The author hereby grants to the U.S. Government and MIT permission to reproduce 
and to distribute publicly paper and electronic copies of this thesis document in 
whole or in part in any medium now known or hereafter created. 

Author. 

Department of Mechanical Engineering and 
Sy stem Desigp and Man agement Program 

May 10, 2019 



Certified by. 



Ihemistoklis P. Sapsis 
Associate Professor of Meclfanical and Ocean Engineering 

Thesis Supervisor 


Certified by 



Bryan R. Moser 

Academic Director and Sr. Lecturer, System Design and Management 

i'hesis Supervisor 


Accepted by 


Accepted by 



Nicolas Hadjiconstantinou 
Chairman, Committee on Graduate Students 
















2 



Quantification of Extreme Event Statistics 
in Ship Design 

by 

Uditbhan S. Rathore 


Submitted to the Department of Mechanical Engineering and 
System Design and Management Program 
on May 10, 2019, in partial fulhllment of the 
requirements for the degrees of 
Master of Science in Engineering and Management 

and 

Naval Engineer 


Abstract 

Increased operational demands on Navy vessels extend time at sea and service life, 
making the accurate prediction of catastrophic failures increasingly challenging. The 
high value of these capital assets puts great pressure on designers and decision-makers 
as they work towards preventing such failures while balancing both engineering and 
material cost. The current method for the quantihcation of extreme events is direct 
Monte Carlo simulation supplemented by complex statistical models. When such 
models are not sufficiently bound by physics-based simulation, the noise of statistical 
uncertainty quickly overpowers the response predictions for rare events. 

This thesis builds on previous work at the MIT Stochastic Analysis and Non-linear 
Dynamics (SAND) lab for the quantihcation of extreme events using wave groups. 
By separating the event probability from the physics models, we are able to capture 
rare events in ship motion and loading conditions for a modest computational cost. 
Improvements to the wave groups methodology ensured the slope and amplitude 
of the incident waves rehected the waves encountered in a given wave spectrum. 
The remaining discussion explores the value of a near-real-time risk analysis tools in 
reference to ship design and ship operations, with unique application to Navy and 
commercial vessels. 
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Chapter 1 


Introduction 


A statistical understanding of ocean waves has been fundamental to the design and 
monitoring of marine and naval structures due to the random nature of the seas. More 
specihcally, the ability to understand the risk and frequency of catastrophic events 
have been topics of interest in numerous helds. In a practical sense, ship designers 
write specihcations to deliver a product that will perform in the worst reasonable 
situation. The main problem they face is that these situations are very rare events; 
otherwise failures would shut down the industry. The proportion of catastrophic 
failures to the number of vessels at sea has shrunk in the past half century because of 
rule-based design standards and improved quality assurance practices. A catastrophic 
failure is in the realm of societally acceptable risk and risk probability, but the fiscal 
and operational cost of losing a capital asset like a Navy warship is still great. For 
example, it is accepted that a combatant naval vessel might be shot at during a 
conflict and the expectation is that within a scope of damage, the vessel should 
survive. Similarly, society expects that a cargo ship will deliver its cargo safely during 
a routine transit, but not necessarily while transiting through a hurricane. Therefore, 
heavy weather guidance systems exist; however, establishing what is acceptable risk 
is not a technical matter, nor in question. Risk that was acceptable 70 years ago is 
not necessarily acceptable today. The more appropriate question, and the question 
of this thesis, is whether a more efficient and accurate approach to the quantihcation 
of extreme events in ship design exists. If so, how would it affect ship design, ship 
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operations, and ongoing science? 

The International Maritime Organization (IMO) took first steps in 2009 to shift 
the maritime indnstry from prescriptive design standards to a risk-based system. This 
shift was driven by the conflict between technological advances in mannfactnring and 
compntation with a design envelope that was limited by regnlation. This foray into 
Risk-Based Ship Design (RBSD) is discussed further in Chapter 2 and provides the 
background for this work on the quantihcation of extreme events statistics on ship 
design. Traditionally, quantifying the risk using physics-based models has been a 
resource-intensive process. Designers and decision-makers get around this challenge 
by augmenting shorter physics simulations with statistical models, accepting a sig- 
nihcant amount of uncertainty in the prediction of rare, catastrophic events. This 
study works to develop a method of characterizing the statistics of the rare events at 
a fraction of the computational cost. 

1.1 Contributions 

This study improves upon Wave Groups methodology developed at the MIT Stochas¬ 
tic Analysis and Non- Linear Dynamics (SAND) lab in an effort towards the quantih¬ 
cation of extreme event statistics of ship motion and loading conditions. A blended 
approach of physics and statistical models was necessary to capture the probability 
and severity of rare events. 

1.1.1 Risk-Based Ship Design 

Scientihc and technological advances and an overall improved technical capability 
have fueled innovation in the shipping sector to meet the demand for larger, more 
specialized, and more capable ships. For decades computational tools have allowed 
designers to replicate model testing, ship systems performance, and model ship pro¬ 
cess in a digital space. By integrating these decision support tools into the design 
process, it is now feasible for designers and decision makers to evaluate a much wider 
design space than ever before. The wave groups approach presented in this work 
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provides a deeper understanding of the probability of extreme events in ship motion 
and loading conditions to inform functional requirements of for a risk-based design 
approach. 

1.1.2 Quantification of Risk in Ship Motion and Loading 

The work presented in this thesis builds directly off previous work completed 
by members of the MIT SAND labs. Sapsis, Stevens and Mohamad provided the 
ground work for the application of the Wave Groups methodology and Probability- 
Decomposition-Synthesis (PDS) that enabled the worked presented here. The em¬ 
phasis of this study was to capture the additional properties of the wave systems that 
produce the extreme events and reproduce them within the wave group simulations. 

1.1.3 Industry Application 

The hnal phase of the study is a discussion on the application of this approach 
for quantifying event statistics using on systems engineering principles. We consid¬ 
ered the technical feasibility and barriers of entry for achieving a near-real-time risk 
decision support system based on wave groups. In the past decades, many attempts 
have been made to measure incident wave characteristics at sea; however, all attempts 
have been overwhelmed by the non-linearities of the sea and the ship’s structure. Fi¬ 
nally, we considered the applications of the Wave Groups methodology for the future 
of ship routing via decision support systems (DSS) driven by operational efficiency 
rather than operating limits, as they are today. 


1.2 Organization 

This thesis is divided in to three sections, using the systems development life 
cycle V-model as its framework. The hrst section, beginning with Ghapter 2, is a 
discussion on risk-based ship design and roll the Wave Groups methodology plays 
in the quantifying risk for ship design and ship operations. Ghapter 6 begins the 
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second section of the discussion, which represents the bulk of the thesis work. Here 
we discuss approach to the analysis and the steps taken to rehne and validate the 
model. Chapter 4 provides details about the simulations setup and the design tools 
used for the analysis and Chapter 5 is the discussion on the simulation and analysis 
is presented in Chapter 6. The hnal section, in Chapter 7, is a discussion on the 
industry application of the Wave Groups methodology. This is presented in Chapter 
7. Concluding remarks and considerations for future work are provided in Chapter 8. 
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Chapter 2 


Risk Based Ship Design 


In this section we discuss the development of risk-based standards in the mar¬ 
itime industry and highlight tools and methods designers have for risk-based decision 
making. Today, advances in computing tools have allowed designers to develop more 
detailed statistical and physics-based models to inform design decisions. Regulating 
bodies have traditionally lagged to keep up with advances in technology and knowl¬ 
edge. We dehne Risk as the product of event severity and event probability. This 
created two major areas of research quantihcation and prediction of extreme events. 
This is the foundation of the work presented in this study. 

2.1 Historical Background 

Most risk analysis approaches can be traced back to the nuclear industry, for 
which risk analysis became important in the 1960s. The goal was to calculate the 
probability of the release of nuclear contaminants. Decades later, the Piper Alpha 
disaster in 1988 forced the offshore industry into risk analysis as a requirement started 
with the Norwegian government and by the United Kingdom in 1992. Many more 
members of the Organization for Economic Co-operation and Development (OECD) 
soon followed [1]. 

Historically, in the maritime industry, regulations were a reaction to major acci¬ 
dents or disasters. In many cases, these regulations resulted from an ad hoc safety 
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assessment process seeking to reduce public and political pressures for action. Such 
processes favor quick solutions over a rigorous technical analysis that explores the 
cost versus benefits of solutions. Worse, this reactive approach has made the regula¬ 
tion system more complex as amendments are continually made, thus leaving limited 
room for innovation and advancement in the industry. Further analysis showed that 
this approach to safety standard was a symptom of prescriptive regulation, where the 
burdened of technical analysis laid on the regulating body and classification societies 
(e.g. IMO and the American Bureau of Shipping), leaving it up to ship builders to 
design within the prescribed constraints. 

The transition to risk-based approaches began with probabilistic damage stability 
standards in the 1960s, but took decades to be introduced by the International Con¬ 
vention for the Safety of Life at Sea (SOLAS). SOLAS is responsible for setting the 
minimum safety standards for the construction and operations of merchant vessels. 
The International Maritime Organization (IMO) maintains and amends these stan¬ 
dards, which, until 2009, were updated as a result to collision accidents and disasters. 
Project HARDER (1999 to 2003) investigated elements of the traditional approach 
and proposed a new formulation for the probabilistic damage stability employing en¬ 
hanced computational models and statistical methods [6]. The final recommendations 
from this study were adopted in 2005, prompting a four-year study. Project SAFE- 
DOR (Design, Operation and Regulation for Safety), across a consortium of fifty-two 
European organizations to incorporate risk-based approaches into ship design, oper¬ 
ations, and regulations [1]. This effort concluded in 2009 when the regulations were 
entered into force at the IMO. 


2.2 Risk-Based Ship Design 

The maritime industry moved to risk analysis or risk-knowledge models as a proac¬ 
tive approach to safety standards to get ahead of future disasters and to improve 
innovation. The maritime industry, in one way or another, has been caught between 
two main design drivers; (1) economic drivers to push more cargo faster and cheaper; 
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and (2) societal drivers to reduce losses at sea and to do less harm to the environ¬ 
ment. Under these pressures and limitations of traditional prescriptive regulations, 
ship designs quickly reached a sub-optimum optimal. Figure 2-1 illustrates the gap 
between prescriptive regulations and technical feasibility [1]. 





'echnicdl feasibility 


Best Solution 


Regulations 

Envelope 


Acceptable 
risk level 
(rational) 


Innovation 

potential 


Mostly 
Unknown 
risk level 


Figure 2-1. Possible design envelopes [1] 


Risk based ship design offers two approaches to realizing the innovation potential 
that is the gap; safety equivalence, where the standard is based on a reference vessel 
or defined by a specified risk acceptance criterion [1]. Both approaches derive from 
the introduction of safety as an objective in the design process. 


2.2.1 Risk-Based Design in the U.S. Navy 

The United States Navy manages some of the world’s most expensive portfolios 
of projects. The average cost of just one of the sixty-six active Arleigh Burke-class 
guided missile destroyers is upwards of $1.8B plus further cost for operation and 
support throughout the life of the ship. The process for the design and acquisition 
of projects of this scale are outlined in the Operations of the Defense Acquisition 
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System, more commonly known as the DoD 5000. Figure 2-2 depicts this Defense 
Acquisition Management Framework. The dotted red lines highlight where in the 
process the hndings of this (and related studies) could be applied using risk-based 
models. 


Phases 


\Atorl< 
Efforts 

Activities 

Technical 



/ 

A 


A9k IOC Foc : 

Concep 

Refiner 

kConcepf 
r Deotior 

t 

lent 

Techno 

Develoi 

ogy 

ment 

System 

Demons 

System 

Integrati 

)evelopment & 
ration 

System 

>n Demonstration 

A Design 
^ Readness 
~ Review 

Production & 
Deployment 

FPIP PrediiUlon 

lOTE aFRP 

AOecMion 

” Review 

Operation) & 
Support 1 
Sustainmeijt 

Disposal 


FTe Systems Acqunition 




Systems Acquisition 
TRA 

▼ 


<— Sustainment » 


Reviews 


A 

ITR 


•A 

ASR 


A — A-A A A 


SRR 

moce 
than (Hte) 

^ Techr>ol<>gy Readiness Assessment 


IBR SFRPDRCDR 
Miybt 
more 

than one^ 

Technical Reviev^ 


A *A — A 

TRR SVRlRytV OTRR 
PRR 


A 

PCA 


ISR 


""A 


A Program Reviews 


Figure 2-2. Defense Acquisition Management Framework [2] 


When designing a ship or writing design specihcation and operating limits, de¬ 
signers and project managers must consider the worst reasonable situation. The main 
problem with these "worst" design situations is that they are very rare and difficult 
to quantify. Risk-based ship design offers a framework to understand the risk being 
accepted and the characteristics of these rare but catastrophic events. The complete 
process requires a combination of physics-based simulations, statistical methods, and 
model testing. During an interview at Naval Surface Warfare Center - Carderock 
Division (NSWCCD), engineers working in the maneuvering and seakeeping division 
agreed with our belief that no combination of these tools is mutually exclusive, rather 
a systematic application of all three should be the recommended approach. 


2.2.2 Model Testing 

Model testing has been a critical aspect of naval architecture since before the 
20th century. Reynolds and Fronde offer similitude principles that allow engineers 
to translate the physics of full-scale operations to a smaller scale and for much lower 
cost. The International Towing Tank Conference (ITTC), established in 1932, is 
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a non-government organization that has been charged with establishing standard 
practices for ship model work, procedures, and methods of result presentations. Today 
their aims include improving methods for, and recommending procedures for, physical 
model experiments, numerical modeling, full-scale measurement of ships and marine 
installations [7]. 

Data collected over the past century of model testing and accidents has created a 
vast wealth of knowledge, which has been used as the basis of computational models 
that allow designers to optimize about a set of design requirements. Mathematical 
models are traditionally used to £11 the knowledge gaps for which testing was not 
or could not be accomplished, and to extrapolate to the far reaches of the design 
spaces that have yet to be explored. Though more economical than full-scale testing, 
model testing is still a slow, resource-exhaustive process that is ill suited for capturing 
statistics of rare events efficiently. At NSWCCD, Ken Weems shared that even after 
three hours of simulation, the mathematical fit may not be right. For instance, a 
Weibull fit for an extreme bending moment model may vary between a factor of two, 
or five in certain cases, dependent on the length of the simulation [8]. When operating 
on the tails, where extreme events reside, conditions change quickly and uncertainties 
have a large impact on the reliability of these models. The belief in this study, is 
that these models can be made more reliable by incorporating more deterministic 
simulations at the extreme cases that lead ship instability. 

Naval architects face two main questions with model testing; (1) how to best 
replicate the real-world wave environment and (2), what are the most interesting 
conditions that are worth replicating? NSWCCD has gone to great lengths to address 
these questions directly. Originally built in 1962, the Maneuvering and Seakeeping 
(MASK) basin allowed engineers to test the performance of ships and moored systems 
in realistic sea conditions using pneumatic actuators. Figure 2-3 is a design drawing 
of the 86.400 ft^ basin. In 2013, the facility completed a six-year renovation project to 
replace the original actuators with 216 individually controlled electro-mechanical wave 
boards. This allowed for greater capability and control of the wave generation process 
at high frequencies, and flexibility to customize the spectral shape and wave-front of 
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the wave systems. As a departure from the traditional tow-tank testing approach, 
very few facilities like this exist around the world that can evaluate maneuverability, 
stability, and control of models. 



Figure 2-3. Historical Drawing of NSWCCD Maneuvering and Seakeeping Basin [3] 

However, to fully characterize the extreme events that a vessel would encounter 
during its twenty-hve-year service life, it would require years of model simulation and 
even longer to make a statistically relevant argument. Computer simulations offer 
a complementary solution to model testing, where by the simulation can be used 
to identify the few seaway characteristics and operating scenarios for engineers to 
replicate and study using model testing. 

2.2.3 Computer Simulations 

Computers simulation offer several alternatives to traditional model testing. Quan- 
tihcation of ship response statistics through transfer function analysis works well for 
the case of linear responses. Solving these problems in the frequency domain, the ship 
response calculations can be resolved with little computational effort [9]. These tools 
have helped the ship design industry to understand the core of the probability density 
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function (PDF). Though useful, this is not necessarily the limiting design case. The 
extreme response phenomenon lay in the tails of these PDFs. The non-linearities 
arise from interaction between the wave and hull geometry. Cousins and Sapsis, at 
MIT SAND lab, have shown that the frequency domain analysis does not sufficiently 
characterize the tail structure of ship roll response [10]. Instead, they noted that 
the extreme the roll events occur more commonly than predicted. Computational 
fluid dynamics (CFD) using potential flow models can reconcile these non-linear in¬ 
teractions and capture these response phenomena. Today, CFD tools substitute a 
signihcant portion of model testing, and in some cases, replace model testing alto¬ 
gether. Potential flow and strip theory provide resistance, stability, motions, and 
loading conditions [9]. To collect data on the short-term statistics for these param¬ 
eters, naval and marine engineers have long used direct Monte Carlo simulation to 
reproduce realistic seaways in the computational environment. Monte Carlo is a 
robust approach, but has a major drawback; the accuracy for capturing long term 
statistics is directly proportional to length of the simulation. Since the input waves 
are selected directly from a wave power spectrum, the probability of a response event 
is coupled directly with the wave spectrum. On a standard PC, the pace of Monte 
Carlo simulation is just shy of being about real-time. This is too slow when the goal 
is to write a design specihcation or standards for an event that might occur a few 
times during the 25 year or longer service life of a vessel. 

2.2.4 Quantification of Extreme Events 

Since the IMO adopted RBSD as a regulatory framework, the maritime industry, 
including the US Navy, has gradually shifted to better quantify and predict rare and 
catastrophic events. The challenge with this approach, lies in modeling these systems 
of systems. In many cases, the physics we use are computationally demanding and the 
cost at stake continues to grow. Consider the Navy’s risk related standards, where 
the requirement is that the ship should be designed for to survive and event that 
happens during 0.0001% of the life of a ship. With an average thirty-year service life 
of a US Navy warship, this equals approximately 200 to 1000 hours of operation [11]. 
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The probability of a catastrophic event outside of this area is beyond the scope of 
the design. To observe such an event, the direct Monte Carlo approach would have 
to produce 500,000 hours of simulation. Further, to achieve the level of statistical 
signihcance, this event should be observed no less than 10 times, totaling 5,000,000 
hours of simulations for a single combination of speed a heading. Quickly, the size of 
this analysis becomes insurmountable even for the most advanced computer clusters. 

The alternative to capturing statistics of these extreme events is what we have 
mentioned already, extreme value theory (EVT) £t by mathematical models. The 
precept for EVT is that rare events are similar in their behavior based on some non¬ 
linear physical property. The second approach is the Wave Groups methodology. A 
wave group is a sequence of waves related by a function of amplitude, period, and 
length scale. The Wave Group methodology has been used in the past to deconstruct 
wave systems into discrete portions. We apply this methodology by capturing the 
deterministic ship response to each wave groups in a large set. By overlaying the 
probability that a ship will encounter the wave group, we are able to calculate the 
statistics of the response event itself. In practice, the EVT and Wave Groups methods 
are complimentary because they seek to solve the same problem using different sets of 
assumptions, for the case of our sponsors at NSWGGD, if both approaches can provide 
the same result, this makes the approval process for a future risk-based standard or 
model more palatable. 

2.2.5 Model-Based Engineering in Ship Design 

The purpose of running these computational simulations is to capture the short¬ 
term and long-term statistics of ship response parameters. This is part of the larger 
vision to quantify the frequency and severity of events to help understand risk for the 
design and decision-making process. Gomputational models offer endless application, 
from understanding complexity to streamlining the design process. Models also pose 
disadvantages. Relying on a model without understanding considering the uncertain¬ 
ties or limitations, can be just as bad as using a random number generator. In all 
cases these system engineering models seek to combine aspects of the physical state 
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with event probability and decisions to approximate a future system state and risks 
associated with those states. 

Project SAFEDOR assembled design teams to tackle the following eight design 
scenarios using risk-based design approaches. In all scenarios, the teams developed 
several innovative tools using model-based system engineering to show how risk-based 
ship design could be incorporated into the design, decision making, operation, and 
regulation on the maritime industry. The impetus of the new design framework 
was for teams to balance performance and cost with a new objective - minimize 
risk [12]. Table 2.1 lists some of the unique applications of risk-based design that 
were accomplished during Project SAFEDOR. 


1 

Optimize a post-PANAMAX size cruise vessel based on new probabilistic 

damage requirements. 

2 

Design a cruise liner that will be safer than existing ones and reduces incidental 

damage to environment in case of grounding or collision. 

3 

Increase cargo capacity of an existing fast full-displacement RoPax ferry. 

4 

Design a RoPax ferry for 50 passengers using risk-based principles vice SOLAS 

requirements. 

5 

Develop a fire risk model to provide a quantitative measure of risk of new designs 

and economic benefit for a new lightweight composite superstructure. 

6 

Develop a short-sea LNG vessel to distribute gas using the safety equivalence of 

the CNG rules. 

7 

Develop a competitive, low-gross tonnage open-top container vessel using safety 

equivalence to overcome the current regulatory disadvantages. 

8 

Optimize a post-PANAMAX size cruise vessel based on new probabilistic damage 

requirements. 

9 

Multi-objective optimization of an AFRAMAX tanker by genetic algorithm to 

improve transport economy and reduce environmental impact. 


Table 2.1. Applications of risk-based design explored during Project SAFEDOR 


As the hrst large scale project to develop the risk-based regulatory framework 
for the maritime industry, the work provided hrst principle approaches to safety and 
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corresponding design tools. During the early stage design, tools for risk analysis could 
be simpler (e.g. expert judgments, databases, or simplified formulas). As projects 
progress in to the concept and detailed design phases, these tools must become more 
rigorous and verihable [1]. Tools such as model-based system engineering and agent- 
based modeling have proven to be able to quantify safety as a measurable objective. 

2.3 Wave Groups 

This study applies the concept of wave groups for the quantihcation of risk in ship 
operation and ship design. The basis of this formulation, also known as critical wave 
groups, is that probability of a certain type of ship instability may be determined by 
the probability of encountering the wave group that generated that instability. This 
approach disassembles the problem in to a deterministic and a probabilistic part [13]. 
This decomposition-recomposition of the problem shown in hgure 2-4 is intended to 
deal with the non-linearity. 



Figure 2-4. Wave group process using decomposition-recomposition approach 


Wave groups have been the focus of much of the work motivated by risk-based 
ship design and risk-related standards. Previous work completed by Anastopoulos 
et al, identihed wave group characteristics such as height, period, and length scale 
through the deterministic evaluation of ship dynamics. A Markov chain model was 
used to capture probabilistic part by determining the frequency of encountering any 
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wave group above a certain threshold [13]. The principal focus of this study was 
to construct wind-generated, high probability wave groups and used them for the 
assessment of ship stability in irregular seas. 

The pitfall is that quantihcation of the wave groups must be consistent with the 
physical mechanisms of the response. Previous work has shown that if we propose 
an arbitrary non-linear response that was only based upon wave elevation, the ap¬ 
proach may work well to capture some of the response parameters, but is unable to 
capture others. Essentially, the problem is more complicated and the wave slope or 
the distribution of the waves over space must be consistent to capture all response 
parameters. 


2.4 A Regulatory Perspective 

The IMO is responsible for publishing stability standards for all vessels. These 
standards provide a general framework for designers. With the focus on safety, these 
guidelines often fall short of maximizing performance and accounting for extreme 
events. 

The rapid pace of technology has brought on design tools and manufacturing 
methods with which traditional prescriptive standards have struggled to adapt. The 
2008 Intact Stability Code published by IMO was the hrst acknowledgement from a 
standards or classihcation body that design technology is rapidly evolving and codes 
should be re-evaluated and revised as necessary. The code further recognized many 
of the challenges with the quantihcation of extreme events that have been discussed, 
specihcally that "complex hydrodynamic phenomena have not yet been fully investi¬ 
gated and understood" [14]. This was a signihcant shift from the retroactive approach 
to a model-based approach. 

Chapter 5.1.1 of the IS Code explicitly states that "compliance with the sta¬ 
bility criterion does not ensure immunity against capsizing, regardless of circum¬ 
stances" [14]. Though understandable, this is a bold statement from the organization 
responsible for setting standards for stability and safety at sea. The problem of quan- 
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tifying the risk due to extreme events at sea remains a challenge that many around 
the globe are actively working to address. The approaches and methods discussed in 
this study represents a novel approach towards this problem. 


32 



Chapter 3 


Theoretical Analysis 


In principal, we know that the formulation of the wave system does not a priori 
consider the formulation of the ship response. Potential flow simulation was used to 
capture the response motion and loading conditions deterministically for both the 
control and wave group experiments. The control data was developed using direct 
Monte Carlo simulation during previous work completed by Stevens [15]. During 
his work, Stevens characterized an original set of wave groups by modulating peak 
amplitude and length scale. Over the course of this study, we sought to characterize 
the slope of the incident wave group. A brief overview of the direct quantihcation 
and wave group implementation is discussed followed by two numerical approaches 
that were implemented for the experimental set. 

3.1 Direct Quantification 

Direct Monte Carlo simulation is one method for reconstructing wave helds from 
a wave power-spectrum. Wave-power spectra differ by regions of the world and are 
developed through extensive data collection over several years of observation and 
experimentation. To implement Monte Carlo, the wave spectrum can be divided into 
equally n-integer subdivisions. The simulation randomly selects the average frequency 
over the subdivision limits based on the probability of occurrence. The wave elevation 
of the resulting linear waves can be superimposed over space and time to produce a 
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realistic seaway using 3.1. An example of the resulting 1-D wave system is shown in 
figure 3-1. 
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n=l 


-1-^^- 

-1-1- 

1 

h't 




1 {1 S,-*) K i", \ A AJ 1 ! j M 

\ ;M GM ' Mi MMi li (m / M 

' ' ■ ■ , ■ Ww V \ 

_^_ 

1 l\k \A'!| V Till'll 

'Mi" li ' ''; 

_1_1_ 

\/ y I'J 1 j ' 1 1 1 ' [1/ 1 / 1/ 1/1 1 W l! \r 

'ij ' ■ ~ 

_1_i_1_i_ 


500 550 600 650 700 750 800 

Time [s] 

Figure 3-1. Wave system developed using Monte Carlo 


A peak wave group amplitude f3 wave heading 
LX peak frequency 9 phase angle 

L length of the wave group k wave number 

In equation 3.1, C, is total wave elevation sum of the linear waves generated by the 
wave-power spectrum. It is generally accepted that the wave elevation of a seaway 
is a stationary ergodic, Gaussian random process [16]. The ocean waves we are 
most interested to reproduce are wind generated wave described by the JONS WAP 
spectrum. More details on the spectrum used for the simulation are provided in 4. 
In summary, the JONSWAP spectrum is narrow banded, as such the wave elevation 
follows a normal distribution and the probability distribution of the peaks can be 
defined by a Rayleigh distribution, 3.2. 
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(3.2) 


3.2 Wave Group Decomposition 

The Wave Groups methodology considers the analytical decomposition of a ran¬ 
dom wave field into short, discrete wave groups. The assumption of a hyperbolic- 
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secant-shaped profile was selected following work completed by Sapsis and Consins 
[10]. The resnltant windowing fnnction r] = A{t)sech{^), shown in hgnre 3-2, was 
implemented to prodnce the modihed wave elevation eqnation, 3.3, [15]. 



r]{x, y, t) = Asech{ — ) cos{k{x cos P + y sin (3) — ut + 6) (3.3) 

Jj 

In eqnation 3.3, we consider a single wave system, whose elevation is Ti{t). The 
wave nnmber is dehned by the dispersion relationship for deep water waves, eqnation 
3.4, where g is gravity, and A is the wavelength. 


k = 


271 

T 



(3.4) 


The key advantage of the wave gronps, is that it can be completely determined 
by a few parameters and the realistic range of these parameters is determined by the 
record we collected from the long Monte Carlo simnlation. Fignre 3-3 is the resnlting 
wave system with the windowing fnnction in red, shown along the peaks of the wave. 


3.2.1 Wave Group by Length Scale 

The hrst implementation of wave gronp sampling was developed nsing the joint 
probability of length scale and amplitnde [15]. This worked well to captnre some 
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Figure 3-3. Incident wave system generated by modified LAMP 


response parameters, minima, maxima and minima, but fell short in cases 
of pitch maxima and minima, F^ maxima, and My maxima and minima. Further 
inspection suggested the quantification of these wave groups was not consistent with 
the physical mechanisms of the response and that the problem was more complicated. 

wave groups reconstructed in figure 3-4 show the wave elevations of two such 
sampled wave groups. Here the peak frequency is constant and the peak amplitudes 
were nearly identical. However, the slope of the wave was not properly characterized 
{blue = 2.878 and red = 3.191). 



Figure 3-4. Wave Elevation Comparison 


Figure 3-5 shows the resulting surge force, Fj,, response. The peak response is 
dramatically different between these two wave groups. The peak force for the 90m 
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wave group is nearly five times larger than the 160m wave group for the same ampli¬ 
tude. The slope of the wave is a key driver in the formulation of the ship response 
and was not a controlled parameter from this first sampling set. 
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Figure 3-5. Wave Elevation Comparison 


3.3 Capturing the Slope 

From the example shown above, we knew that the wave slope needed to be an 
input in the form of probability distribution rather than a product of the windowing 
function. However, we were limited in the inputs that we could give LAMP to L, A, 
Uo, X, and </>. To extract the frequency from the time series data, we considered a 
simplified, 1-D, wave system with a random phase angle, 0 = [0 : 27r]. This simplified 
the incident elevation equation for the wave groups to 3.5. 



(3.5) 


The gradient of the rj with respect to time is the slope of the time series 3.6. 



(3.6) 


If we consider the slope and elevation to be independent variables, we can use 
their product to develop a joint distribution map. From this we can calculate the 
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frequency, cu, for individual wave groups using the approach in 3.7. 


rj' max[ujAsech{^)sm.{kx — ut)] 
T] max[Asech{^) cos{kx — ut)] 


(3,7) 


The hyperbolic secant and amplitudes canceled, and the maximum value of one for 
both sine and cosine, also canceled. The remainder was the frequency with units of 
radians per second. 



V 


(3,8) 


The CO- variance of the one million data point set of slope and amplitude was 0.85. 
However, after considering multiple smaller sets, we found that the co-variance failed 
to converge to a single value, in other words, what we were observed was just noise. 
This was enough to conclude that slope and elevation were indeed independent sets. 
This was made simpler because both the wave elevation and its derivative followed 
Rayleigh distributions, as expected. 
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After normalizing each distribution to one, we can calculate the joint probability 
density function as the product of the two distributions. 


fh,h'{0i,92) — fh'{9i) ■ fh{92) 


(3.11) 
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Chapter 4 


Experiment Setup 


This chapter describes the major components of the experiment, to include the 
simulation software, the hull geometry, and the input seaway characteristics. 

4.1 Simulation Tool 

4.1.1 Large Amplitude Motions Program (LAMP) 

This study primarily used the Large Amplitude Motions Program (LAMP) to 
simulate the wave-body interaction problem. Developed in 1989 by the Science Appli¬ 
cations International Corporation (SAIC) and later, the Leidos Corporation, LAMP 
calculates time-domain motions and loads for floating bodies using 3-D potential-flow 
panel method and viscous effects are included as an external, non-pressure forces in 
the overall force calculation. The challenge with any physics-based formulation that 
incorporates nonlinear wave-body interaction is in balancing the computation require¬ 
ments with correctness and complexity. To do this, LAMP offers the user multiple 
options for the complexity of calculations. 

• LAMP-1; Body-linearized 3-D 

• LAMP-2; Approximates large-amplitude 3-D body-nonlinear method with as¬ 
suming steady forward speed 
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• LAMP-3: Approximates large-amplitude 3-D body-nonlinear method allowing 
large lateral motions 

• LAMP-4: Large-amplitude 3-D body nonlinear method 

The LAMP-2 method sufficiently captured the body nonlinear effects for this study 
with slightly more computer resources than LAMP-1. This method used the body- 
linear 3-D approach to compute the hydrodynamic part of the pressure forces, while 
the nonlinear hydrostatic restoring and Froude-Krylov wave forces were calculated on 
the actual hull surface below the incident wave surface. The very important aspect of 
the time-domain approach implemented in LAMP is that it provides and event-driven 
analysis of the ship response to a specific wave event rather than a general statistical 
solution of the ship motion problem. 

4.1.2 LAMP System 



Figure 4-1. Structure of LAMP-Based Analysis System 
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The LAMP System is a comprehensive suite of software tools that include pre¬ 
processing, hydrodynamic calculations, general data processing, flow visualization, 
and finite-element structural analysis. Figure 4-1 from the LAMP manual explains 
the structure of a LAMP-based analysis. 

This study used LAMP for direct Monte Carlo simulation as well as the wave 
groups sampling experiments. Previous work by Stevens developed MATLAB code 
that automated the process of generating input files, executing LAMP simulations, 
and parsing the output data for statistical analysis [15]. This study improved upon 
this process by providing more tools to stream line the process while implementing 
additional methods of characterizing wave groups. Figure 4-2 illustrates the imple¬ 
mentation of the LAMP system and MATLAB for this work. 



MATLAB 

Modified LAMP code 


Figure 4-2. Implementation of LAMP and MATLAB for this study. 


4.1.3 Modified LAMP 

The original LAMP source code as adapted to accommodate the Wave Groups 
methodology discussed in 3. An explanation of the modification is provided in 
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Appendix A. This new version, which has come to be known as Modihed LAMP 
(mLAMP), was developed Stevens while at MIT SAND lab [15]. Modified LAMP 
allowed us to completely define the wave group with three variables, frequency, am¬ 
plitude, and length scale. 

4.2 Input hull characteristics 

One advantage of LAMP’s 3-D physics-based formulation is that can accept both 
conventional and multihull structures. The formulation, however, is intended for dis¬ 
placement hulls and is not suitable for planing or surface effect ships. In coordination 
with on-going work MIT SAND Lab, we selected the ONR tumblehome hull for which 
the geometry and wave groups data that was readily available. The model was de¬ 
tailed enough to have research value, but was not representative of any U.S. Navy 
hull in service, allowing the results of this study to be released to the public. 

The ONR tumblehome hull is representative of ocean-going naval surface combat¬ 
ant vessels. The tumblehome characteristic was most recently seen in the US Navy’s 
Zumwalt-cl&ss destroyer. Advantages of this hull include superior radar signature re¬ 
duction, as well as improved acoustics due to the reduced wake from the wave piercing 
bow. A primary disadvantage of this hull, and the reason it has been a topic of study, 
is the decreasing waterplane area due to the inward sloping freeboard. Table 4.1 lists 
the specific hull characteristics for this hull. 


Characteristics 

Value 

LBP 

154 m 

Beam 

18.78m 

Draft 

5.49 m 

Depth 

14.5 m 

Displacement 

85071 

Block Coefficient 

0.535 


Table 4.1. ONR Tumblehome hull characteristics 

lamp’s graphics-based utility program, LMPLOT, was used to generate three- 
dimensional renderings of the input hull. Figure 4-3 shows the ONR tumblehome hull 
that was used for this study. 
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Figure 4-3. ONR tumblehome hull views 



Figure 4-4. ONR tumblehome isometric 3-D view 


4.3 Input seaway characteristics 

Wind-generated ocean wave systems are best characterized by their power-spectral 
density functions. LAMP can reproduce seven input seaways: no seaway, multiple 
linear waves superimposed, a single non-linear Stokes wave, a long-crested 2nd order 
irregular Stokes wave, the Brettschneider spectrum, the Pierson-Moskowitz spectrum, 
the JONS WAP spectrum and the Ochi spectrum [17]. 

It was appropriate to continue to use the Joint North Sea Wave Observation 
Project (JONSWAP) ocean-wave spectrum to validate the new Wave Groups method- 
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ology against previous work [18] [15]. The JONSWAP spectrum was originally se¬ 
lected for its well-studied, irregular ocean spectrum. It accounts for non-linear wave- 
to-wave interactions that make it a suitable spectrum of long-term statistical analy¬ 
sis [16]. The LAMP expression for the JONSWAP spectral density is: 


S^if) 


ag 

/-y5 


-5 exp 


^ / '^m \ 4 


(4.1) 


ojm = 27r/p the peak frequency 

5 = 0.07/or/ < fp 

6 = 0.09/or/ > fp 

The LAMP input parameters for JONSWAP are: fetch ( 7 ), signihcant wave 
height(iLs), and modal period ( 0 ;^)- To effectively compare the results of this study 
with work completed by Stevens, the JONSWAP parameters remained unchanged 
between the two works. The resulting spectral density plot is displayed in 4-5. 



Omega (rad/s) 

Figure 4-5. LAMP input JONSWAP spectrum 

JONSWAP Spectrum: a = 0.060,7 = 3.0,77^ = 13.2 m, and Tp = 10 sec. 

The long time-series wave elevation closely followed a Rayleigh distribution as 
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expected 4-6 




Wave Elevation Maxima 



Figure 4-6. Wave Elevation Statistics 


4.3.1 Data Corrections 

A MATLAB script was used to count the wave peaks that followed mean up- 
crossings and the maximum slope in of the peak. These peaks are marked by a 
red circle and were used to develop the wave elevation maxima PDF. Nearly 8% of 
peaks occurred between — 1 m to 1 m. In these cases, it was difficult to reconcile the 
incident frequency with the approach discussed in Chapter 3, resulting in frequencies 
that were abnormally large relative to the realistic waves systems. Moving forward, 
we considered only peaks whose elevation was greater than Im and noticed little 
change to the wave elevation PDF. 

To obtain the 1000-hour time series, 2000 30-minute LAMP simulations were 
concatenated. This was necessary to avoid self-repeating wave systems which is a 
known limitation of LAMP [15]. Each LAMP simulation also contained start-up 
and stop transients that skew the ships motions and loads. The start-up transient 
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ramps all wave elevations between factors of 0.0 and 1.0 over the first 120 simnlation 
time-steps. When processing the data, the first 300 time-steps and the last 100 time- 
steps were removed to ensnre accnrate results. However, this method introduced 
discontinuities into the time-series data that needed to be accounted for in wave 
slope calculations, figure 4-7. 



Figure 4-7. Discontinuity in time series data 


Reproducing the time series data without discontinuities was possible, but is out¬ 
side the scope of this work. The alternate solution was to exclude all slopes that 
exceeded 13.8, which is 0.003% (1099 points out of 3.5 million). The slope cutoff of 
13.8 was chosen because it coincided with the point where Monte Carlo PDF diverged 
from the Rayleigh distribution, figure 4-8. 



Wave Slope [dh/dt] 

Figure 4-8. Monte Carlo wave slope PDF 

The Rayleigh distribution verified that the peak slopes were representative of a 
physical wave system since we know the derivative of a Gaussian process, should also 
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be Gaussian. 


4.3.2 Generating Probability Maps 

We started with the probability map of the previously discussed JONS WAP spec¬ 
trum that Stevens used to produce the wave groups of various amplitudes and length 
scales, hgure 4-9 (left). Even though we shifted to a new approach, it was still neces¬ 
sary to establish an input length scale for the wave group. Figure 4-9 (right) shows 
the surge force response for these wave groups [15]. We noticed that force was rel¬ 
atively constant for 6*i greater than 2 or 45 m. The value of = 0.0402, from the 
10 sec peak period. 



[Lengthscale / k^] b H-engtuscaie i 


Figure 4-9. (Left) Length-Amplitude Probability Map. (right) Surge Maxima Response Map 

This phenomenon occurred in all response cases which allowed us to choose a 
single length scale value of 3.6 or 90 m. This length scale was long enough to reduce 
variability in the results and short enough to keep simulation run-time as low as 
possible. By selecting a single length scale, we were able to maintain the two degrees 
of freedom analysis and maintain the simplicity of the problem. 

Figure 4-10 shows the hnal distribution of the slope of the peaks and the joint 
probability map of the slope and peak amplitude. The probability map was combined 
with the wave groups sampling simulations to produce response maps and response 
distributions. 
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Figure 4-10. (Left) Wave slope {drj/dt) PDF. (right) Slope-Amplitude Probability Map 


4.4 MATLAB Analysis 

The LAMP systems offered an extensive suite of tools to generate and evaluate 
large sets of simulations. However, given the duration of this study, MATLAB pro¬ 
vided considerable capability with a much smaller learning curve. Also, this work was 
built on previous work completed by Stevens and Mohamad who had also developed 
many tools in MATLAB for pre- and post-processing LAMP simulations. 

4.4.1 Wave Groups Pre- and Post-Processing 

As previously mentioned, the work completed by Stevens [15] provided data for 
1000 hours of simulation. He used MATLAB to prepare, execute, and concatenate 
the short-time series results to produce this data set. From this set we were able to 
generate the probability maps using the procedures discussed in Chapter 3 for the 
wave groups analysis. 

This study developed adapted versions of these tools to generate distributions for 
the wave slopes and frequencies that fed into the wave group simulation. A major 
portion of this work also included optimizing and streamlining these tools to provide 
greater functionality and be robust to user error. 
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Figure 4-11. Progress plot output 

4.4.2 Wave Simulation and Analysis Tools 

During the development and testing of the simulation and analysis methods, 
it became quickly apparent that validating a large set of wave groups would be a 
monotonous ordeal. The sampling sets grow as the square of the sampling density. 
As we sought high resolution, the sampling density grew to 25, and later 50, restarting 
a 74-hour simulation after 50 hours was not sustainable. Additionally, as part of the 
quality assurance and understanding of the simulations it was necessary to hlter and 
select interesting cases to animate. 

The hrst tool develop was a simple script, ProgressPlot.m to read the hies in the 
LAMP output directories to check for completed simulations. The script provides the 
user with a plot of completed wave groups by folder, 4-11. With this knowledge, a 
new set simulation could be started where the previous one hnished (allowing gaps to 
be readily hlled along the way). Furthermore, having the hexibility of knowing where 
each simulation resided allowed us to divide the simulations over multiple directories 
preventing corruption or total loss of a large set. 

We found a simple solution in the analysis method. Since we were already using 
the wave groups to generate a response map for motions and loading conditions of 
interest, we could use resulting response maps as an interface to select the points 
that most stood out from the sampling set. The ResponseMap.m function gave the 
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Figure 4-12. Selection of a single sampling simulation using the ResponseMap.m function. 


LAMP file Z_p_o_s_A14.68_dh0.8761_f0.05968.out 



0.00 190.0 0.00 190.0 

T= 134.60 

Figure 4-13. LMPLOT simulation of selected wave group sample 


user the option to select a point on the response map of interested, which would be 
copied into a separate directory for further study. A white circle indicates the selected 
simulation and the red crosses indicated lines of constant frequency. Figure 4-12 is 
an example of an application of this click function. 

The screen capture of the resulting simulation shows the failure of this LAMP 
simulation. This result was subsequently removed from the sampling space, Figure 
4-13. 
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4.5 Response PDF Computation 


Once the wave group sampling set was completed, we manually consolidated the 
output hies into a single directory. The LAMPOut.m script removed the start-stop 
transients from each response time series, hgure 4-14 and selected the maximum and 
minimum value, hgure 4-15. 



Figure 4-14. Processing LAMP transients 



Figure 4-15. Wave Group max & min 


These values were stored in an n-by-n matrix to produce simulation response maps. 
Once each simulation was processed, the script combined the wave group probability 
map and the response map to produce the response distributions. Chapter 5 discusses 
the simulations performance in more detail. 
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Chapter 5 


Wave Group Simulations 

5.1 Computational Tools 

Fast computational performance has been a prerequisite to achieve practical appli¬ 
cation of simulating a hydrodynamic environment. Much of the work at MIT SAND 
lab has been to hnd new and innovative ways of reducing this computational cost. 
Although computational time depends largely the machine being used, the times 
recorded provide a representative view of the current performance achieved in an 
academic environment. All LAMP simulations and data processing was completed 
on a single desktop PC, Intel Core 17-3770 at 3.40 GHz, 16 GB RAM. 


5.2 Simulation Performance 

Early simulations attempted were run with as few wave groups across the prob¬ 
ability map as possible. The since size of the wave group sampling set grows by 
the square of the number of samples across each distribution. The lower resolution 
sets, proved the simulations ran, but did not adequately produce reliable distribu¬ 
tion curves of the ship response and loading conditions. Table 5.1 shows how the 
computational time changed with simulation size. 
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Sampling Density 

3 

5 

10 

15 

25 

50 

Number of Wave groups 

9 

25 

100 

225 

625 

2500 

LAMP Simulation Time [hrs.] 

0.27 

0.75 

3.0 

6.7 

18.6 

74.3 

Data Processing Time [min] 

0.21 

1.0 

2.4 

5.3 

14.8 

59.8 


Table 5.1. Simulation Summary 


By evaluating the full range of amplitudes and slope distributions, we were able 
to update probability maps independent from the modihed LAMP simulations. The 
LAMP simulations made up the 96% of the 78 hours of the computational time. The 
time spent reprocessing response for a new probability map was a fraction of the 
total time, table 5.1. The sampling size of 50 provided enough resolution to produce 
a smooth distribution for the response analysis. Previous work by Stevens showed 
that a greater sampling density provided little to no improvement in results [15]. 


5.3 Simulation Outputs 

5.3.1 Motion and Loading Response 

Figure 5-1 is a short time series set of all the response parameters from the original 
Monte Carlo simulation. By simulating wave groups that approached directly off the 
bow, we were able to simplify the analysis from 12 response parameters to hve, table 
5.2. All values were calculated about the ship’s center of gravity (CG). Evaluating 
motions caused by beam or quartering seas was outside the scope of this study and 
has been left for future work. 
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Figure 5-1. Time series summary of response parameters 


1 

Vertical position of CG [Zpos] 

2 

Surge force at CG [Fr] 

3 

Heave Force at CG [F^] 

4 

Pitch about CG [Yrot] 

4 

Pitching Moment about CG [My] 


Table 5.2. Response Parameters of Interest 


5.3.2 Failed Simulations 

Figure 5-2 is the run time for each slope-amplitude wave group simulation. By 
tracking simulation time, we learned that run time was a good proxy for simulation 
performance. The average simulation processing time was 107 seconds, so simulations 
that ran for less than 70 seconds, ended early because of a capsize or simulation failure. 


55 

















































Simulation Run Timefs] 
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Slope [m/s] 


Figure 5-2. Slope-Amplitude simulation run time 


Figure 5-3 shows the two types of simulation failure observed during the analysis. 
The capsize, left, is the result of an extremely low slope and large amplitudes. The 
plot on the bottom right corner shows a 14.5 m drop in wave elevation, causing the 
ship to rapidly lose stability by losing waterplane area and thus, capsizing. The failure 
to the right shows a vessel capsizing because of encountering multiple 14 m to 16 m 
waves at a frequency of 12 s. 



Ta M7.00 


Figure 5-3. Types of simulation failure. 
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5.3.3 Response Map Adjustments 


Though the intent was to reproduce realistic ocean waves, we found that certain 
combinations of amplitude and slope produced non-physical wave groups in LAMP, 
resulting in extreme response calculations that made it difficult to analyze the remain¬ 
ing portions of the response maps. Figure 5-4 shows the original response map for 
the heave force minima (left) and adjusted response map (right). The color gradient 
shows the change in the minimum measured heave force for each wave group, slope, 
and amplitude combination. 



'2 4 6 8 10 12 Xio® 

Slope [m/s] 


Minima [kN] 



2 4 6 8 10 12 

Slope [m/s] 


Figure 5-4. Improvement in by filtering out transient spikes. 


The overlaid on the response map are contours the joint wave group probability, 
hgure 4-10. Notice the outlier at the far left of the unadjusted heave force response 
map. Further inspection of this simulation showed an unexpectedly large start up 
transients, which were not related to the wave groups themselves, figure 5-5. 
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LAMP -File F_x_Fl0.60‘1Z_dh5.382_f8.907.out 



r-^ 

tal(l) 




0.00 Eta 320.0 


-280 



Figure 5-5. Unexpected heave response at low wave slope 


Similar cases were present in other response parameters. To improve the resolution 
in the areas of interest, (center and top right of the response map), we smoothed 
out the effects of transients for slopes less than 0.7m/s. This did not impact the 
distribution calculations, but allowed us to see the rest of the response map with 
greater detail. Further discussion on the physical meaning of these response maps is 
provided in Chapter 6. 
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Chapter 6 


Analysis of Simulation Results 


As briefly discussed in Chapter 5, results from the wave group sampling simula¬ 
tions were collected in response maps. Of the twelve output parameters considered, 
flve produced interesting results for further analysis. In this chapter we discuss each 
of the response parameters individually. 

The wave groups act as a threshold value, but to establish speciflcation for op¬ 
erational limits or risk-based standards, it was necessary to consider the wave group 
probability map. Finally, we compared the PDFs of the slope-amplitude sampling 
simulations with the previously completed lengthscale-amplitude approach, as well 
as the Monte Carlo simulations. The forces, moments, and motion measured in the 
following simulations are with respect to the ship center of gravity. 


6.1 Surge Force 

Amplitude-Slope Wave Groups 

The surge force, Fj., is the forward-aft force felt due to accelerations along the 
X — axis of the ship. The ship was given a steady heading of 10 kn, so additional 
accelerations were introduced except those caused by the wave group. Figure 6-1 is 
the response maps from the slope simulations. The contour bands represent the joint 
probability of the wave group combination in a logarithmic scale. 
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Figure 6-1. Surge Force Response Map of the Slope Simulation 


Bands of constant heave force are clearly visible along the diagonals of the slope- 
aniplitnde response map. These bands align with lines of constant freqnency, which 
were a product of calculations discussed in 3. The two regions where the simulation 
failed or capsized are also present in these maps. The filtering step, reduced the 
impact of the extreme values on the far left of the map, but the region at the top 
of the map with extreme values is still present. This extreme response area has very 
distinct edges resulting from resonance between the hull and the large wave group. 
Figure 6-2 shows a snapshot of the simulation highlighted by the white circle on the 
response map. 


LAMP file F_x_A15.64_dh7.634_f0.4882.out 



21000 

0.0 


0.00 

T= 143.50 


-14.0 

0 



Figure 6-2. Wave Group Simulation A =15 m ^ =7.63 ms ^ 
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This is an extreme sea way, and likely outside the operational limit of most sea go¬ 
ing vessels. This wave is 15.6 m (50 ft) above the mean waterline and is categorized as 
Sea State 9 on the Douglas Scale, by the World Meteorological Organization (WMO). 
As noted by the probability contours, such waves systems are very rare and result 
from extreme hurricane conditions, with sustained winds over 67 kn. For reference, 
the bridge of an Arleigh Burke-class Destroyer sits 55 ft above the waterline. 

LAMP file F_x_All.8_dhl.877_f0.1591.out 



-owu I _ I_ _I -i.u.ui _ k _ i 

0.00 320.0 0.00 320.0 

T= 138.20 


Figure 6-3. Wave Group Simulation A =12 m ^ =1.88 ms ^ 

For comparison, we considered a more modest wave group, shown in figure 6-3. 
This is still a relatively large wave system, with peak amplitude of 12 m (39 ft), but 
the peaks are much further apart than the overall length of the ship, so the extreme 
conditions goes almost unnoticed. 

Response Distribution 

We combined the response and probability maps from the sampling set to develop 
response distribution. The slope-versus-amplitude sampling set was compared to the 
Monte Carlo simulation, show in blue in figure 6-4, and lengthscale-versus-amplitude 
sampling set, shown in yellow. The middle row is in natural scale, while the bottom 
row is provided in a logarithmic scale in the y-axis. The Rayleigh distribution is 
also provided as a dashed line to show the Gaussian approximation. Since we are 
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considering only "extreme events," the PDFs are only provided for cases beyond two 
standard deviations from the mean valne (fnrther referred to as the "tail"). The 
gray vertical lines depict the mean (thick grey) and standard deviations (increasingly 
lighter gray lines). 




Figure 6-4. Surge Force Distribution 


The slope-amplitnde sampling set consistently over predicted the ships response 
in all loading conditions. The Rayleigh distribntion nnder predicted the minimnm 
snrge force, snggesting a heavy tail strnctnre. The lengthscale-amplitnde sampling 
set remained the best fit for both responses. The slope simnlations least accnrate 
followed the shape of the Monte Carlo simnlation. 


6.2 Heave Force 

The heave force, F^, is the vertical force along the z — axis of the ship. Fignre 
6-5 is the response maps from the slope-amplitnde simnlation set. 


62 











































2 4 6 8 10 12 

Slope [m/s] 


Maxima [kN] 

■ i ftW 



2 4 6 8 

Slope [m/s] 



Figure 6-5. Heave Force Response Map of the Slope Simulation 


For this response, we selected a simulation that ht a more realistic design limit 
within the hrst contour line. Highlighted by the white circle in the minimum surge 
force map, 6-6 shows 6.6 m (21ft) wave system. This heave force is on the order 
of 100kN (lit), which would still be categorized as an extreme event for the ship 
structure, though, it did not result in a stability failure. 


LAMP file F_z_A6.682_dh3.129_f0.4683.out 



T= 



6.750 


-6.75 


Eta 



0.00 


320.0 


Figure 6-6. Wave Group Simulation A =6.6m uj =0.47rads ^ 


Moving forward, providing additional snapshots and response maps did not add 
signihcant value to this discussion, but are available for reference in Appendix B. 
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Response Distribution 


Figure 6-7 shows PDFs for the heave force minimum and maximum. The length 
scale (LA) distribution performed best in the minimum force distribution. No specihc 
distribution stood out as htting the Monte Carlo distribution better in the maximum 
force distribution. 



Time [s] 



Force [kN] xio'' 



Figure 6-7. Heave Force Distribution 


6.3 Heave Motion 

Heave motion was measured as the vertical position of the ship center of gravity 
taken from its initial position at f = 0. Figure 6-8 shows the PDFs of the sampling 
sets. These results show no clear correlation between any wave group sampling set 
and the Monte Carlo distribution. 
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Figure 6-8. Heave Motion Distribution 


6.4 Pitch Motion 


Pitch is the rotation about the transverse or y-axis of the ship measured in degrees. 
The upward motion of the bow relative to the center of gravity is considered positive 
pitch or maxima. Figure 6-9 shows the PDFs of the sampling sets. The slope sampling 
more closely captured the shape of the Monte Carlo distribution than in previous 
response cased. In both the maximum and minimum, it was interesting to note that 
the slope sampling and the Rayleigh distribution, over predicted the weight of the 
tails. 
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Figure 6-9. Pitch Motion Distribution 


6.5 Pitching Moment 


The final parameter we considered was the pitching moment, My. The minimnm 
moment is in moment in the negative, bow-down direction. Fignre 6-10 shows the 
PDFs for the wave group sampling sets. The slope-amplitude sampling approach was 
able to capture the shape of the Monte Carlo simulation only in the My minima case. 
Pitch and pitching moment were the main parameters we sought to capture with a 
new wave group sampling approach. 
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Figure 6-10. Pitching Moment Distribution 


6.6 Summary 

The slope-amplitude sampling distributions over predicted the response in all 
cases. Rather than capturing the maximum value of the response, we can bring 
the distribution closer to the Monte Carlo simulation by taking the average of the 
three or four highest values. This would be more representative of the entire wave 
system. Also, we can improve the sampling efficiency by taking into account wave 
steepness. Such high frequency wave systems do exist in nature but generally exotic 
waves. We had expected the joint distribution maps to make the effects of these 
combinations negligible, but ultimately, these systems impacted the final response 
parameters more than predicted. 
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Chapter 7 


Industry Application 


The aim of this analysis was characterizing the extreme events in terms of a series 
of short simulations that may be used for further analysis using high hdelity CFD or 
basin model testing. Though several discrepancies remained in our characterization 
of wave groups, a successful realization of this approach would provide a basis to 
know exactly the risk of each wave system and physically what the loads were going 
to be [8]. In this way, we would not only be able to reproduce the conditions that 
result in these rare events, but also understand the conditions associated with them. 

Wave groups were incorporated for the hrst time, into the Second-Generation 
Intact Stability Criteria. This shifted the basis of stability criteria from accident 
statistics to physics-based. The new criteria are formulated in tiers, with the lower 
tier being focused on early stage designs and the follow-on stages requiring additional 
evaluation of "safety level" i.e. the probability of failure despite having satished the 
criterion [13]. 

7.1 Routing for Operational Efficiency 

The successful application of the wave groups concept could provide ship operators 
and owners with a risk evaluation tool that incorporates operational loads into navi¬ 
gational routing decisions. Similar Decision Support Systems (DSS) already account 
for fuel efficiency and heavy weather avoidance. Figure 7-1 illustrates the decision 


69 




Figure 7-1. Route optimization for fuel efficiency 



Figure 7-2. Route optimization for heavy weather avoidance 


model for route planning on the basis of fuel efficiency. The iterative loop considers 
route waypoints and the operating prohles of the main engines to determine best 
course and speed. This is generally the hrst phase of route planning. 

Concurrent to fuel routing, national weather services like the Nation Oceanic and 
Atmospheric Administration and the Navy Weather Center, provide weather update 
using national weather assets and observational data provided by other mariners. 
This information is generally passed along as sea state recommendation that can be 
correlated to the vessehs operational limits. These limits are predetermined dur¬ 
ing £rst-of-class rough water trials or through response amplitude operators (RAO) 
simulations. Figure 7-2 illustrates this decision model. 

This is useful for pre-voyage route planning and providing updates along the 
way. However, this approach is limited to avoiding extreme excursions outside the 
operational profile. Papanikolaou et. al argued that for on board DSS to be truly 
practical for route optimization and risk mitigation, probabilistic assessment methods 
would need to replace traditional deterministic guidance models [19]. The Wave 
Group methodology offers just this. Weather data can be translated into wave group 
characteristics and overlaid on the navigational chart. Using a library of a pre¬ 
computed, standard set of wave groups and associated responses, the DSS is able 
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Figure 7-3. Route optimization for operational efficiency 


to forecast operational loads directly on the navigational chart. Figure 7-3 illustrates 
this decision process. Ultimately, the operator sees a list of potential routes and 
filters. The filters are based on operational criteria that can reduce life-time fatigue 
and improve engine efficiency. 

7.1.1 Stakeholder Analysis 

We identified two primary stakeholders that would benefit from a such a DSS: 
The Navy and the commercial maritime industry. However, we found that their 
needs and how they would implement this system would be very different. Consider 
the operating profiles of a standard Naval combatant vessel and a loaded bulk cargo 
vessel, figure 7-4 [4,5]. Navy ships operations are based on a larger mission set 
of training and deployments. During a deployment a vessel may need to rapidly 
reposition from one area to another and then loiter in that area until operational 
commanders provide further tasking. 

In the case of commercial vessels, the operational tempo, though still demanding, 
is much more predictable. As shown in this speed profile, once the vessel has left 
port, the route and speed are well established. Commercial vessel operators also have 
the option to delay shipment in case of heavy weather or re-route to reduce extreme 
loads. For these reasons, both stakeholders are interested in very different ends of the 


71 



























40% 


40% 


i. 35% 

30% 
S 25% 

S 20% 

E 

'■S 15% 
o 

C 10% 

o 

^ 5% 
^ 0 % 


















Navy Combatan 

t 
















L998 

Yearly Speed 1 

Profill 

e 




























































































/ A \^ \ 





/ 

















1—^ 


0 5 10 15 20 25 30 

Speed (kts) 



Figure 7-4. Annual speed profile for a US Navy combatant (left) and a loaded bulk cargo vessel 
(right) [4,5] 


trade space, figure 7-5, and realize the benefits of the Wave Group methodology and 
near-real-time risk evaluation in different ways. 
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Figure 7-5. Operational efficiency routing tradespace and stakeholders 


Navy designers seek to better understand the characteristics of extreme events so 
that they may account for the demanding operational loads from the conception of the 
design. Wave groups offers a physics-based model to replace existing mathematical 
models and reduce uncertainty for the extremely rare end of the distribution. In the 
end, if the current estimates over predict certain loading conditions, the specihcation 
can reduce plate thickness and weight. Alternatively, in the case current estimates 
under predict extreme events, actions can be taken to reduce further damage and 
potentially mitigate a total loss of a $1.8B surface combatant. 
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The commercial maritime industry, is a much large stakeholder in terms of shear 
numbers of vessels and annual miles traveled. The Wave Groups methodology opens 
and new type of routing for operational efficiency, which is applicable even in the 
absence of heavy weather. By correlating loading conditions to incident wave patterns, 
ship operators can develop routes that avoid the heaviest loads and reduce fatigue 
for their specihc hull structures. Consider the transit between two of the busiest 
ports in the world, the Port of Shanghai, and the Port of Los Angeles, which are 
approximately 5,633 nautical miles apart [20]. Though modeling this fatigue analysis 
was beyond the scope of this work, the benehts could possibly add years to the service 
life of these vessels. In 2018, over 1800 vessel arrived at the port of Los Angeles, of 
which 73% were from China and neighboring countries [21]. The value of such an 
approach is only multiplied by these factors. 

7.1.2 Technology Risk 

With current technology. Navy and cargo vessels could be sufficiently instrumented 
to track loads and fatigue imposed on the hull during normal operation. This data 
would then be used to improve the current wave groups theory and vice versa, a 
successful wave groups model could identify deviation in normal operational loads as 
means of predicting failure. The estimation of sea state parameters on the moving 
ship remains a fundamental problem to be solved [19]. Wave groups also offer an 
alternative to full scale wave prediction by reducing the scope to the prediction of 
smaller sets; however, this work is still theoretical. 

Nearly a decade ago, attempts to capture wave elevation in real time used X-band 
radars to map the wave elevation and while simultaneously recording the non-linear 
motion of the ship with accelerometers. Unfortunately, the ship does not perform like 
a buoy and many efforts like this over the past decade have tried and failed to reconcile 
this problem. The most recent developments of sea state sensing were completed 
from a shore facility in Istituto per 1‘mbiente Marino Costiero (lAMC) of the Italian 
National Research Council in 2017 [22]. Other recent efforts have applied machine 
learning and neural networks to help quantify extreme events in ship response, but 
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training these systems has proved challenging because the ocean environment changes 
too drastically from the training area to the operational area. In the interim, the 
National Oceanic and Atmospheric Administration (NOAA) continues the densely 
monitor the world‘s oceans. Wave elevation by buoy and wind data are still the best 
predictors of heavy weather. 

7.2 Risk Mitigation with Digital Twin 

An innovative application is in the development of a shore based digital twin. 
If the wave group approach can quantify risk in near-real-time, operators and ship 
owners have the opportunity, via a digital twin, to determine the real risk of a vessel 
at sea. This would provide a more dehnitive way of understanding risk beyond sea 
states and polar plots, which are used today. Wave groups would provide a cata¬ 
log of precomputed responses to the wave systems and in real time understand the 
conditions the ship is in and what is the probability of response events. 

7.3 Ship Design and Model Testing 

Once an approach is accepted, reproducing the wave groups in the Maneuvering 
and Seakeeping (MASK) Basin is straight forward. Instead of running weeks of model 
tests in the basin to evaluate a new design, the design team can run the set of 
waves known to lead to instabilities [11]. With the probability of these waves groups 
occurring in a certain sea state in a given area, these sets can be easily accommodated 
into heave weather guidance systems. 
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Chapter 8 


Conclusions & Recommendations for 
Future Work 


An approach to the construction of wave groups by characterizing slope and ampli¬ 
tude are presented in this study. In all cases, the response distributions over predicted 
the Monte Carlo simulation, which could be caused by several aspects data process¬ 
ing approach. It is possible that rather than selecting the largest response, using the 
average of a few of the most pronounced peaks could bring down the tails. From the 
sampling perspective, a more rigorous method for filtering out non-physical waves 
using parameters such as steepness could be implemented prior to the mLAMP sim¬ 
ulations. This would reduce the number of simulations required, saving memory and 
time. Several MATLAB tools were provided to aid in the continued analysis of wave 
groups and mLAMP simulations. The final results of this wave group study were 
compared to direct Monte Carlo quantihcation and an amplitude-length-scale char¬ 
acterization of wave groups. Neither slope nor lengthscale approaches stood out as 
exclusive solutions to characterize the non-linear responses. The vertical position re¬ 
sults in both cases were completely inconsistent, it would be worth exploring how this 
data was collected. Overall, the Wave Group methodology not only shows promise 
for both the quantification of extreme events, it also has potential the potential to 
improve tools for operational efficiency routing. 
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8.1 Future Work 


Several opportunities for future work exist. The slope-amplitude probability map, 
consistently over predicted the Monte Carlo distribution. Rather than tracking the 
peak response amplitude, it would be worth taking the average of the top three or so 
peaks for each wave group time series rather than just the maximum value. There is 
an implementation of a mean crossing function that was used during the full Monte 
Carlo data analysis that would need to be implemented in the LAMPout.m script 
for each response parameter. 

This and the previous application of wave groups to quantify ship response statis¬ 
tics simplihed the probability space into two dimensions. Stevens characterized the 
wave groups by lengthscale and amplitude. Here, we characterized the wave groups 
by slope and amplitude. A hnal attempt would re-capture the lengthscale distribution 
to produce a three-dimensional probability space. This a heavier problem computa¬ 
tionally, but this effort could reduce that burden with a more deliberate sampling 
approach. Instead of evenly spacing samples across the PDF. It is possible to spread 
out samples at low probability regions of the PDF, while also hltering out non-physical 
or exotic wave groups. 

One of the applications, in discussion with the project sponsors, was the improve¬ 
ment in model testing time. Future work might consider combining the LAMP wave 
group simulation with model testing at the MIT tow tank. 
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Appendix A 


LAMP Modification 


Pefromed by Stevens in 2017 and provided as a reference [15]. 

One of the LAMP seaway inputs is for a superposition of sinusoidal waves with 
surface elevations of the form [17]; 

N 

r]{x,y,t) = An COs{kn{x cos/3n + ysm/3n) - (Xj- + 9n). 

n=\ 

These waves are similar to the desired form, and were identihed as a good starting 
point. The primary modihcation was the inclusion of the windowing function S 
displayed in Figure 


V*{x,y,t) = S{x) *C{x,y,t)- 


S{x) sech{^) exp(^)+exp(f)‘ 

The author obtained the LAMP source code and identihed the hies governing the 
wave input functions. The windowing function was added into the incident wave 
elevation function. While the experiment calls for a single wave input, the function 
was left in the summation form to avoid unnecessarily modifying the code. The 
incident wave elevation was therefore modihed to: 
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N 

r]{x,y,t) = AnSech{-) COs{kn{x cos/3n + y sin f3n) -Unt + 9n). 

n=l 

As long as only a single incident wave input is recieved, N will be equal 1, and 
this equation would revert to the desired input function; 

X 

ri{x, y, t) = Asech{—) cos{k{x cos [3 + ysin P) — ut + 9). 

Jj 

The length scale input into the modihed version of LAMP was separate from the 
wave input hie, therefore only a single length scale can be used even if the analysis of 
several linear waves is deemed important. 

As stated previously, LAMP does most of its calculations using the potential 
function. Within the superposition of sinusoidal waves seaway input, the potential 
function is calculated as: 

^{x,y,zA) = ^^exp''"^sin(A;„(a;cos/3„+ |/sin/3„) -Unt + 9n). 

n=l 

This function also needed modihcation to include the windowing function; 


^*{x,y,zA) = S^{x,y,zA)- 

Which for a single sinusoidal wave input can be simplihed to: 


<h*(a;, y, z, t) = —sech{—) exp^^ sin(A;(a; cos /3 + y sin/3) — ut + 9). 

CO Jj 


Additional modihcations were made to the time and spatial derivatives of the 
potential function 0 and incident wave elevation function rj. The result of these 
modihcations was 
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Appendix B 


Complete List of Simulation Figures 
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Figure B-1. Surge Force Response Map - Slope 
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Figure B-2. Surge Force Distribution 
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Figure B-3. Heave Force Response Map - Slope 
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Figure B-4. Heave Force Distribution 
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Figure B-5. Pitching Moment Response Map - Slope 
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Figure B-6. Pitching Moment Distribution 
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Figure B-7. Pitch Response Map - Slope 
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Figure B-8. Pitch Distribution 
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Figure B-9. Heave Motion Response Map - Slope 
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Figure B-10. Heave Motion Distribution 
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